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Properties of the charmonium-like state X(3872) are investigated and its nature 
ON ' is discussed as based on the existing experimental data. In particular, we analyse 
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the new data from Belle and BaBar Collaborations and argue that, while the BaBar 
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data prefer the dynamically generated virtual state in the DD* system, the new 
Belle data clearly indicate a sizable cc 2 3 P\ component in the X wave function. 

PACS numbers: 14.40. Gx, 13.25.Gv, 12.39.Mk, 12.39.Pn 



I. INTRODUCTION 

In a few recent years a number of new states in the spectrum of charmonium have been 
found experimentally. These states, labelled as X's, K's, and Z's attract special attention 
of phenomenologists since most of them (if not all) can hardly fit into the standard quark 

o 

On 
O 



model scheme. This means that, in addition to the genuine cc component, the wave functions 
of these states must have extra components, whose nature is not yet clear and is an open 
problem. Although various scenario are suggested and discussed in the literature, such as 
threshold phenomena, hadronic molecule, and so on, no unambiguous criteria which would 
allow one to distinguish between different assignments for these "homeless" charmonia have 
been established so far. 

Among these new charmonium-like states the X(3872) meson is most well-studied. The X 
was first observed in 2003 by the Belle Collaboration in charged S-meson decays B — > KX, 
in the mode 7T + 7T~J/ip with the dipion originated from the p-meson. The mass and the 
width of the X(3872) reported in Ref. [l[ were 

M x = 3872.0 ± 0.6(stat) ± 0.5(syst)MeV (1) 

and 

r x < 2.3MeV. (2) 
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Later Belle reported their observation of the same state in the tt+tt 7i J/ip (uJ/ip) and 
jJ/ip modes [2], with branching fractions 



Br(B -> KX)Br{X -> 7J/VO = (1.8 ± 0.6 ± 0.1) x 10~ 6 , 
Br(X -> 7r + 7T-7r°J/^) 



Sr(X -> 7T + 7T- J/V») 

5r(X -> 7 J/V») 



1.0 ±0.4 ±0.3, 
= 0.14 ±0.05. 



(3) 
(4) 

(5) 



Br(X -> tt+tt- J/V) 

The X(3872) was confirmed in the discovery mode by the CDF |3], D0 |4J, and BaBar 
Collaborations. In their recent updates, BaBar [6] and Belle |z] reduced slightly the 



branching ratio, 



Br(B + -> X+X) 5r(X -> tt+tt" J/tjj) = (7 ± 10) x 10" 



(6) 



and both Belle and BaBar also observed the X in the 5° decays with the rate comparable 
with the charged channel 6|, |7| • In addition, BaBar measured other decay modes of the X 



Br(X -> -fj/jj) 
Br(X ->• tt+tt- J/V) 

5r(X -> 7^') 



0.33 ±0.12, 



5r(X -> TT + TT" J/lp) 

and imposed the upper limit on the X production [9|: 



1.1 ±0.4, 



(7) 
(8) 



Br(B -> XX) < 3.2 • 10" 4 . (9) 
The most recent CDF result for the mass of the X observed in the tt+tt-J/^ mode [h]] 



is 



M x = 3871.61 ± 0.16 ± 0.19 MeV. 



(10) 



Finally, the quantum numbers J PC = 1++ are favoured for the X (although 2 h are not 
yet excluded) [11]. 

Clearly the measured properties of the X(3872) raise a number questions concerning 
its nature. In particular, although the X is produced in the 5-mesons decays with the 

^that is with the branching ratio typical for 



branching ratio of order 10~ 4 — see Eq. ([9 
genuine charmonia (such as J/ip, ip', or Xci) 
of a 3 Pi cc meson in the vicinity of the observed mass of 3872 MeV (see, for example, 13]) 



12J, quark models fail to predict the existence 
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In addition, quark-antiquark interpretation of the X faces a further challenge, namely, a 
strong isospin violation — see Eq. (j3J). 

In the meantime, the observed mass of the X(3872) measured in the 7r + 7r~J/?/> channel, 
appears to be quite close to the position of the D°D*° threshold which, according to the 
most recent CLEO data 14j lies at 



M D o D * = 3871.81 ± 0.36 MeV. 



(11) 



It is quite natural to assume then that the X wave function contains a large admixture 
of the DD* molecule component 1 , and the isospin violation is readily explained as due to 
the large (about 7 MeV) mass difference between the charged and neutral DD* thresholds. 
Indeed, if one assumes that the decays X — > pJ/ip and X — > uJ/ip proceed via DD* loops, 
then the isospin violation happens due to the difference between the charged and neutral 
loops which, in turn, is due to the aforementioned mass difference. Although it is not large 
per se, it is enhanced due to kinematical reasons, as the effective phase space available in 
case of the p is much larger than that in case of the u 15j, [l6j. 

To summarise, data seem to indicate a dynamical origin of the X. From the theoretical 
point of view several assignments for the latter are discussed in the literature. 

It was noticed long ago that one-pion exchange can be responsible for the formation 
of near-threshold states in .D-meson systems. In particular, one-pion exchange is attractive 



in the 1 
Ref. 



DD* channel 



0, Q, Q, 



201 ] . The X(3872) as a virtual state was discussed in 



16j, as generated dynamica 



On the other hand, in Ref. 



ly from the interaction of pseudoscalar and vector mesons. 
2l| , in the framework of a coupled-channel microscopic quark 
model with the cc — DD* mixing, the X(3872) is generated as a virtual state in the DD* 
channel together with the 2 3 Pi charmonium resonance. Similar phenomenon is found in 
a coupled-channel analysis [22|, where a more sophisticated QCD-motivated approach to 
light-quark pair creation is developed. 



The X as a loosely DD* bound state is advocated in Refs. |23|,[24(, while the cusp scenario 
for the ir + ir~J/ip excitation curve in the X(3872) mass range is discussed in Ref. 25]. 

Meanwhile, although the molecule assignment for the X seems to be quite plausible, it 
meets certain obstacles as well. To begin with, a natural worry is that in the DD* system, 



1 An obvious shorthand notation is used here and in what follows: DD* = -^=(DD* + DD* 
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bound by the one-pion exchange, the pion may go on-shell and thus binding may not be 



strong enough 



151 ] . For the most recent work on the possibility for one-pion exchange to 



discussed in Refs. 



bind the DD* system see 26|, [271]. Further implications of the nearby pion threshold are 



28 



29]. 



Furthermore, quite a large branching ratio for the radiative decay X — > jip 1 (see Eq. (jS 
can be explained naturally in the framework of quark models. Indeed, it is well-known 
that in so-called Coulomb+linear quark potential models the radiative decay x'ci ~^ l^Jjk 
is suppressed in comparison with decay x' c i ~^ 1$' • For example, the estimates of Refs. [13] 
and [30|] yield: 



r(x'd -> 7«VV0 = 70 keV [13|, 11 keV [3 



r(x'd 7^0 = 180 keV [13(, 64 keV [3 



(12) 



In the molecular model an opposite pattern was found in Ref . [19| , and a large 'ytp' rate is now 
considered as an evidence against the molecule interpretation. Notice, however, that there 
exists a mechanism for the radiative decays of molecules via -meson loops which was not 
considered in Ref. [3] and which favours the x'ci ~^ 70' decay rate over the x'ci ~ > jJ/^ one - 
However, a reliable evaluation of such radiative decays of molecules meets severe problems 
with divergent loop integrals which can hardly be resolved in a model-independent way. 



Finally, for a pure molecule, the branching fraction B — > KX was estimated in Ref. 31] 
to be less than 10~ 5 , that is much smaller than the experimental data on the X production 
(though, being very model-dependent, such estimates should be treated with caution). So it 
seems quite reasonable to assume that this is the cc component of the X to be responsible 
for the X production in B meson decays and for the X radiative decays. 

The interest to the X(3872) was catalysed even more in 2006, when Belle reported an 
enhancement of the D°D°7r signal just above the D°D*° threshold observed in the reaction 
B + -> K + D°D 7i° m 



33| at 



M x = 3875.2 ± ± °- 8 MeV 



(13) 



with the branching 



Br(B+ K + D°D 7i°) = (1.02 ± 0.3ljftg) x 10" 4 . 



(14) 
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The corresponding state was called the X(3875) and it was confirmed later by the BaBar 



Collaboration as well 



34]. 



Although an immediate and the most natural conclusion is that this is simply yet another 
manifestation of the same well-established state X(3872), a 3 MeV shift in the mass may 
have had dramatic consequences for such an interpretation. As a result, a rather extreme 
assumption was made that two different charmonium-like states might reside in the same 



mass region. It was noticed in Ref. 



35], however, that, under certain assumptions on the 



nature of the X, the two states could be indeed reconciled with one another. In particular, 
it was argued in Ref. 35j that the X, as a virtual state in the DD* system, can reproduce 
both sets of data, for the DDn and p J j ip channels, whereas in the latter case one deals with 
a threshold cusp. Parameters of the model were tuned to fit all the data on the resonance 



width and branching ratios 
with the findings of Refs. 



21 



'hese results of Ref. 



351 ] appear to be in a good agreement 



251 ] . The analysis of Ref. 35j is improved in Ref. 36| where 



additional non-DD* modes of the X were taken into account and the admixture of the 
genuine cc charmonium in the X wave function was estimated. 



Recently Belle Collaboration announced a new analysis for the D*°D° case 37]. The 



new data on the D 



*o 



D 



°7T 



and D 



*o 



D 0/ ~f channels were fitted both with the simple 



Breit-Wigner line-shape form and with the Flatte formula. As a result, a lower position of 
the peak, 

M x = 3872.6+^;l ± 0.4 MeV, (15) 



was obtained than the one reported before — see Eq. (U3l) . The corresponding branching 
ratio was measured to be 



Br(B + -> K + X(D D )) = (0.73 ± 0.17 ± 0.13) x 10" 



(16) 



These new data are analysed in Ref. 38j using the technique very close to that of Ref. 35] , 
and the conclusion is made that the A is a 2 3 Pi cc state strongly distorted by couple channel 
effects. 

In this paper we present an updated Flatte analysis, with the new data on the _D°D* 
mode 37] and on the jip' mode jsj included. In particular, we address the question of a 
possible x'd charmonium admixture in the wave function of the A(3872). The strategy 
employed in this paper, differs significantly from the one of Ref. 35] , where a model-blind 
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Flatte analysis was performed. Here we assume a mechanism for the X production via the 
charmonium component. 

The paper is organised as follows. In Section [III we give necessary details of the 
Flatte parametrisation for a near-threshold resonance and apply this technique to the case 
of the X(3872). We analyse the data in Section HTT1 and comment on the effect of the D*° 
finite width in Section 0VJ We conclude and discuss the results in Section [V] 



II. FLATTE PARAMETRISATION 



In this Section we introduce a Flatte-like parametrisation of the near-threshold ob- 
servables related to the X(3872) state. Let us define the energy E relative to the neu- 
tral D°D*° threshold (see Eq. ffTTj) ). Then the relevant energy range is approximately 
-10 MeV < E < 10 MeV which covers both the three-body D°D 7i° threshold at 
E D o DOn o « -7 MeV and the charged D+D*- threshold at E D+D *- = 5 « 7.6 MeV. A 



natura 
nance 



generalisation of the standard Flatte parametrisation for the near-threshold reso- 



39j of the D°D*° scattering amplitude reads 

1 9ih 



F(E) 



with 



D(E) 



E — Ef 
E — Ej 
E — Ef 



2k 1 D(E)' 
g X K X g 2 n 2 , ,r(E) 



2 



+ 1- 

1 1 — — h 



9ih 



2 

g2k 2 



+ 



and 



E < 
< E < 5 
E>5 



(17) 



(18) 



Ki 



y/-2mE, k 2 = ^2fx 2 (E-S), k 2 = - E). 



Here \L\ and fi 2 are the reduced masses in the D°D*° and D + D*~ channels, respectively. 
Assuming isospin conservation we set g\ = g 2 = g. 

The term ir(E)/2 in Eq. (|T8l) accounts for non-DD* modes: 



r(E) = r 7T + 7T -j/ 1 p(E) + r 7T + 7T ~ n oj/^(E) + r , 



(19) 
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where we single out the first two modes because of their explicit energy dependence: 

-M-mj^ dm q{m)r p 



r, 



7T+7T - J/lfl 



(E) = S P 



r, 



7r+7r _ 7r° J/lf) 



(E) = k 



2vr ( m _ mp )2 + r p V4' 

M ~ m W dm q(m)r u 



3m x 



2tt (m-mJ 2 + r 2 /4' 



(20) 
(21) 



with f p and being effective couplings and 



q{m) 



(M 2 - (m + m J/v ,) 2 )(M 2 - (m - m^) ; 
4M 2 



(22) 



being the centre-of-mass dipion/tripion momentum (M = + M(D°D* )). 

Now, if we assume the short-ranged dynamics of the weak B —> K transition to be 
absorbed into the coefficient B, then the differential rates of interest in the Flatte approxi- 
mation read: 



o n*o\ 



and 



dBr(B KD°D 
d~E 

dBr(B -> Ktt+tt-J/iIj) 
-> Kir+7r-ir°J/ip) 



-- B 
B 



2tt l^^)! 2 ' 

(^) 
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2tt 1^(^)12 ' 

7r _ 7T° J ftp (E) 



(23) 
(24) 

(25) 



dE 2tt I^D^)! 2 

Obviously, the rate ( |23l) is defined for > only, while the rates (1241) and ( |25l) are 
defined both above and below the D°D*° threshold. Strictly speaking, one is to take into 
account a finite width of the D*, which is very small however. Indeed, the total width of 



the D* ± -meson is measured to be 96 ± 22 keV 



12]. There are no data on the D*° width, 



but one can estimate the total width of the D*° from the data [12j on charged D* ± to be 
about 63 keV, which gives r(D 



*o 



D 



42 keV. If, nevertheless, the finite width of 



the D* is taken into account, the rate ( J23l) continues to the region E < and interference 
effects are possible in the final state, as described in Ref . [40] . We shall discuss this in some 
detail below though, at the moment, we follow Refs. 35|, |36( and neglect the D* width. 
Now, in order to proceed to the branching ratio to the DDir final state, one is to take 



into account the branching fractions of the D 



*o 



12 



Br(D*° -> /yV°) = (61.9 ± 2.9)%, 



Br(D*° -> £>°7) = (38.1 ± 2.9)%, 



(26) 
(27) 
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so that 

dBr{B -> KD°D%°) r ^ 1 g/c 



dB a628 ^W (28) 

Analogously we have for the D°D 0/ y differential rate: 

^ = °- 38S ^^yF' ( } 

With the Flatte parametrisation introduced above, one can make use of the method 



suggested in Ref. 



4l| to estimate the admixture of a bare x'd state in the wave function 



of the X. Indeed, in the context of cc-DD* coupled-channel model the quantities entering 
the Flatte-type expressions for differential rates acquire clear physical meaning. Namely, 
the coefficient B can be viewed as the branching fraction B — ► Kx' c i, g is the bare x'c\DD* 



coup 



ing constant, and i~o is the bare total width of the x' c \ level- Moreover, as shown in 



Ref. 41], in the Flatte limit, the probability w(E) to find the bare state in the wave function 
of a physical state can be expressed in terms of Flatte parameters as 

W{E) = 2^\D{E)\^ 9klQ{E) + 9hQ{E ~ 6) + F{E)) - (30) 

The admixture W of the x'& charmonium in the resonance wave function can be defined as 

/■s max 

W= w(E)dE } (31) 

J -Emm 

where the integral it taken over the near-threshold region. As was discussed before, we 
choose is to be from -10 MeV to +10 MeV. 



III. DATA ANALYSIS 
A. Essentials and constraints 



In this chapter we analyse the existing data using the Flatte approach described above. 
Our aim is to estimate the admixture of the charmonium component of the X wave function 
and to identify the nature of the residual, dynamically generated part of its wave function. 
In particular, we shall answer the question as to whether the existing experimental data are 
compatible with the bound state or virtual state in the D°D*° system. 

Let us comment briefly on the difference between the bound/virtual state situations. 
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In the Flatte approximation the inelastic differential rate is 

b 27r\D(E)\ 2 ' [ ] 

One can see easily that the behaviour of the inelastic rate below threshold depends strongly 
on whether there is a zero in the real part of the denominator D(E) below threshold. Indeed, 
if D(£'bound) = for some inbound < then, for r(E) — * 0, the inelastic rate (l32l) does not 
vanish, but becomes proportional to a ^-function: 

B (^mP~\ E E ) 5 ( E ~ E ^). ( 33 ) 

We end up therefore with a real bound state, which is not coupled to inelastic channels. 

On the contrary, if there is no such zero (virtual state case), the rate (|32|) vanishes as 
r(E) —>■ 0, while the D°D*° rate does not vanish in this limit. 

Consider the case of r = first. In this case, in order to distinguish between these two 
scenarios (bound state versus virtual state) one is, as was argued in [24|, [35(, to check the 
ratio 

Br{X -> tt+tt-J/^)' 

which varies from quite small (and hardly resolvable experimentally against the background) 
values, for the bound state scenario, up to values of order ten (in Ref. 35] this ratio was 
calculated to be 9.9) for the virtual state. It follows from the data quoted in Eqs. flBJ) and 
( fl4|) that the ratio ( |34|) is indeed large (~ 10 4- 15), which seems to indicate the virtual 
state nature of the X. However, the above consideration was based on the assumption that, 
once produced, the X state can only decay through one of the three channels: D D°tt or 
pJ/ip and ioJ/ijj (the 'jJ/ip mode is small, and was neglected). Nowadays, a new jip' mode 
is observed. Moreover, if it is presumably due to cc bare seed, then extra decay channels 
typical for charmonium should exist for the X, which are encoded in the extra width r ^ 0. 
These are annihilation modes (into light hadrons), and Xci(3515)7T7t (the latter was estimated 
in Ref. [42] to be of order of a few keV). The total width of the Xci(3515) is 0.89 ± 0.05 
MeV, and the branching fraction into radiative jJ/ip mode is about 36% [12]. If it were 
a true guide, then one expects the width of the x'& t° be about 1 -j- 2 MeV. Quark model 



prediction 



30] yields the value of 1.72 MeV for the total width of the x' c v ^ n accordance 



with predictions (fl2|) . radiative modes are not the dominant ones. 
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Now the ratio (|34|) should be modified to read: 

Br(X -> D°D°it ) 



(35) 



5r(X -> nornD /) ^) 

opening the possibility for the X being a bound state 2 . 

Finally, assuming the X to be produced via the x'd component of its wave function, one 



can estimate the coefficient B. The world average for the Br(B + — ► i^ + x c i) is 12] 



5r(5 + -> K+Xci) = (5.1 ± 0.5) x 10" 4 , (36) 

and it is known [12] that J/ip and '0' are produced in the B ^ K decays with comparable 
branching fractions: 

Br(B + -> K + J/tp) = (10.22 ± 0.35) x 10 -4 , 
Br(S+ -> K + ijj') = (6.48 ± 0.35) x 10" 4 . 

Then it is reasonable to assume that the x' c i * s produced in the B — ► K decays with the rate 
comparable to (1361) . There exists a quark model prediction 44j Br(B — > -ft'Xci) = 2 x 10~ 4 . 



However, the model used in Ref. |44j underestimates the rate ( 1361) more than two times. 

As was mentioned before, the admixture of the genuine charmonium in the X wave 
function is given by the quantity W defined in Eqs. (1301) and (I3T1) . 

Therefore, our analysis strategy is to approximate the existing experimental data on the 
DD* and 7T + 7T~J/ijj decay modes of the X with the Flatte formulae and 

• to find the admixture of the x'i c charmonium in the X wave function by evaluating 
the integral (|3T!) of the spectral density (|30|) over the near-threshold region; 

• to compute the scattering length for the DD* system and thus to make a conclusion 
concerning its virtual/bound state nature; 

• to investigate the effect of the finite width r . 

The data on the DD* and tt + 7t~ J '/ip modes are analysed under the following constraints: 



2 The idea that, including an extra width, one can fit the data on the X(3872) both with virtual and 
bound state was first presented in Ref. 43 1. 
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. HQ. 



TABLE I: The sets of the Flatte parameters for the Belle data from Refs 

Set T g E f , MeV f p / w B • 10 4 <j) W BW ■ 10 4 a, fm 

1 1.1 0.3 -12.8 0.00770 0.04070 2.7 180° 0.19 0.5 -5.0 - il.3 

2 1.0 0.137 -12.3 0.00047 0.00271 4.3 153° 0.43 1.9 3.5 - il.O 

3 2.0 0.091 -7.8 0.00090 0.00523 3.7 152° 0.52 1.9 3.3 -il. 7 



• Br(B — > Kx'd) = & = (3 -r- 6) • 10 4 , with the preference to lower values (see Eq. ( |36i) 
and the discussion following it); 

• Br(B -> iifX) = iW < 3.2 • 1CT 4 (the limit imposed by the BaBar data js], see 
Eq. ©); 

• i~o = 1 -T- 2 MeV (as discussed above). 

Throughout this paper we deal only with the data on the charged £>-meson decays, as 
the uncertainties in the data on the neutral mode remain large. Belle Collaboration presents 
the data on the D D°ir and D°D 0/ -f modes separately, and we analyse only the former mode 
(again due to larger uncertainties in the D°D 0/ -f mode). BaBar data presented are for all 
D°D*° modes, so we consider these data. 

B. Belle Collaboration data 

As it was mentioned in the introductory part, recently Belle Collaboration presented a 
new analysis for the ir + ir~J/ilj 1 DDn, and DD^y decay modes of the X [3, [s^]. These new 
data differ significantly from the old ones. The peak in the tt + tt~J/iI) mass distribution 



is shifted to the left, making the virtual state/cusp scenario advocated in Ref. [35( less 
plausible. However, as the ratio ( 134"1) remains large, extra non-DDrr modes are needed in 
order to arrive at the bound-state solution, as it follows from Eq. fl35|) and will be shown 
below. 

In order to translate the differential rates into number-of-events distributions, we notice 

that there ar e 13! sig nal events ht the Belle data lot the ^ W > *a„ne. g which 

corresponds to the branching fraction of about 8.1 • 10~ 6 ; the bin size is 2.5 MeV. Then 

131 \ dBr(B -> Kn+n-J/^j) 
8.3 • lO" 6 j dE 



N^<*(E) = 2.5 [MeV] ( ) — ^ ^ " tlZl . (37 ) 
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FIG. 1: Differential rates for the ir + ir~ J/ip channel (left plot) and D°D°-k channel (right plot) 
(see Eqs. ([57]) and (|38|). respectively) with the parameters given by set 1 (see Table H}. The 
distributions integrated over the bins, with resolution function taken into account, are shown as 
- dot8 , .p— ^ (s ee M . QH) are gi ve n - o P e„ dotS w lth error ta 

Similarly, for the D°D n° mode, the Belle Collaboration states to have 48.3 signal events 
in the charged mode 371 ]. which corresponds to the branching fraction of about 0.73 • 10~ 4 ; 
the bin size is 2 MeV. Thus the number-of-events distributions is calculated as 

48.3 \ dBr(B -> KD D°ir ) 



<iif 7r °(^) = 2.0[MeV] 



(38) 



0.73 • 10" 4 / dE 
In the latter case, the background function is proportional to the two-body D°D*° phase 
space i?2 oc yE, that is the background is considered to be due to the contribution of the 
D°D*° and, as such, to interfere with the signal: 



dBr(B -> KD°D°n 
dE 



fci 

0.62— 

2tt 



tie — , + ccosi 



+ Im- 



+ c sin i 



(39) 



D{E) V V D{E) 

with the relative phase 4> and c being fitting constants. 

Finally, the resolution functions for both reactions are taken in the form of Gaussians 
with the fixed resolution scale being 3 MeV, for the n + ir~J/ip channel, and with the variable 
mass-dependent resolution function a(m) = ay/m — mo, with a = 0.172 MeV 1//2 and tuq = 
M(DD*) 0- 



The Belle data on the _D u _D°7r° mode 



371 ] can be equally well described by both the virtual 



state and the bound state in the DD* system (set 1 and sets 2, 3 in Table fl] and plots in 
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FIG. 2: The same as in Fig. [I] but for set 2 (upper plots) and set 3 (lower plots). 

Figs. [1] and [21 respectively). For set 1, the width r mimics the jif)' decay channel and is 
fixed through the condition that Br ('yip') ~ BriimJ/ip). However, the description of the 
Belle data on the ir + n~J/ip mode is remarkably poor for this set — see Fig. [TJ Besides that, 
the radiative width of about 1 MeV seems to be suspiciously large and, in any case, is not 
compatible with x' c i assumption. 

We find therefore that a decent description of the Belle data on the it + it~ J/ij) mode is 
only possible with a bound state (see sets 2 and 3 in Tableland in Fig. [2]). Furthermore, 
because of a considerable contribution of the finite width i~o to the spectral density w(E), 
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its integral over the near-threshold region appears to be rather large (see Table [T|) which 
indicates a significant admixture of the genuine charmonium state in the X wave function. 

It is instructive to study the behaviour of the spectral density for the bound-state case 
in more detail. It is plotted in Fig. [3J where the contribution of non-DZ)* modes is also 
shown which peaks at the position of the bound-state mass. Using the relation between the 



Flatte parameters and the effective range parameters established in Ref. [4l|], it is straight- 
forward to demonstrate that, in the limit of vanishing inelasticity, the spectral density below 
the D°D*° threshold becomes, similarly to the inelastic rate (1321) proportional to a ^-function, 

w{E) -> Z5(E - £ bound ), E < 0, (40) 

with the coefficient Z being nothing but the famous Z-factor which was introduced by 



Weinberg in Ref. |45|] and which defines the probability to find a bare state in the wave 
function of a physical bound state with the binding energy inbound- So it is reasonable to 
define an integral over the near-threshold region: 

Z = / w inel (E)dE, (41) 



En 



with 



Winel(E) = 2n ^ m2 nE). (42) 

Then the factor Z can be viewed as the Z-factor of our bound states smeared due to the 
presence of the inelasticity and it takes the values: 

Z = 0.31 (set 2), Z = 0.37 (set 3). (43) 

The values of the branchings Br(B — > Kx' c i) — B an d Br(B —> KX) = BW, as given 
in Table [U agree with the constrains imposed on them by experimental data and quoted 
in the beginning of this section. The radiative decay width appears to be in a reasonable 
agreement with quark model estimates (J7J): 

r( 7 V>') = 60 keV (set 2), r( 7 V>') = 110 keV (set 3). (44) 

Therefore, one can make the conclusion that the new Belle data favour the X(3872) to be 
a mixture of a genuine charmonium and a dynamically generated molecule-like state which 
appears to be a bound state of the D°D*° system. 
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FIG. 3: Spectral density for set 2 (first plot) and set 3 (second plot). The full spectral density w(E) 
is plotted with the solid line, the contribution Wi ne i{E) of the tvotv-DD* channels to the spectral 
density is shown with the dashed line. The functions w{E) and Wi ne i{E) coincide below threshold. 

C. BaBar Collaboration data 



In the data analysis procedure, similarly to the Belle data case, the formulae for the 
numb er-of-e vents distributions were taken to be: 



^BaBar iP) = 5 [MeV] — - n _ fi — , (45) 



93.4 \ dBr{B -> Kti+ti-J/^) 
8.4- 10- 6 / rf^ 

for the 7T + 7T~J/tp mode (bin size is 5 MeV, number of events is 93.4, and Br(B — > 
Kn+Ti-J/tp) = 8.4 ■ lCr 6 — see Ref. and 

<£» = 2.0[MeV, (46) 



^1.67 ■ io- 4 y 

for the D°D*° mode (bin size is 2 MeV, number of events is 33.1, and 5r(5 — > KD°D*° 



34J). The signal-background interfer- 



1.67- 10 -4 ; all D°D*° modes are included — see Ref. 
ence is taken into account in the same manner as for the Belle data — see Eq. fl39|) . with 
the factor 0.62 omitted. 

The resolution function for the 7t + tt~J/iJj channel is taken in the form of a Gaussian with 
the fixed resolution scale being 4.38 MeV p. As to the DD* resolution, it is described by 
the BaBar Collaboration as a very complicated function, and it is not available in public 
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TABLE II: The sets of the Flatte parameters for the BaBar data from Refs. a, [34]. 



Set 


r 


9 


E f , MeV 


f P 


fu, 


B- 10 4 





w 


BW • 10 4 


a, fm 


4 


1.0 


0.225 


-9.7 


0.0065 


0.0360 


3.9 


113° 


0.24 


1.8 


-4.9 6 


5 


2.0 


0.145 


-6.0 


0.0040 


0.0230 


3.6 


109° 


0.34 


0.8 


-3.9 -i2 


6 


1.0 


0.080 


-8.4 


0.0002 


0.0010 


5.7 


0° 


0.58 


3.3 


2.2 -iO. 3 


7 


2.0 


0.090 


-9.0 


0.0005 


0.0029 


5.5 


0° 


0.53 


2.9 


3.3-i0.7 



domain. In the present analysis we take, with corresponding reservations, this resolution 
also to be Gaussian, with the resolution scale of 1 MeV. 

The BaBar D°D*° data 34j are very similar to the old Belle ones 33], while the ir + 7r~J/il) 
peak in Ref. [f| is moved a bit to the left in comparison with the old BaBar data on the 
same reaction, and the peak width has decreased around 25% due to a better resolution. 
One expects therefore, that the D°D*° data are better described as a virtual state, while 
the 7r + 7r~ J/ip data complies better with the bound state. Correspondingly, we employ two 
different analysis strategies. First, we reconcile the 7r + 7r~ J/ip and D°D*° peaks with each 



other, as it was done in 



35 



(sets 4 and 5). The second strategy is to find the best overall 
description of the both data sets (sets 6 and 7). The parameters for these sets are given in 
Table [TT1 and the differential rates are shown in Figs. H] and [5j 



The reconciling procedure yields a virtual state, similar to the one found in Refs. [35l. 136]. 
The admixture of the charmonium is not large for these sets and, again, the radiative decay 
width r( / -ftp') seems to be too large for the charmonium assignment: 

r(7V>') = 800 keV (set 4), r( 7 ^') = 500 keV (set 5). (47) 

The overall description of the data seems to be not too bad, as the present resolution cannot 
confirm or rule out the 7t + tt~ J/ip cusp scenario. 

Solutions 6 and 7, obtained from the overall fit to the BaBar data, clearly prefer the 
bound-state, similarly to the ones given by sets 2 and 3 for the Belle data. The charmonium 
admixture is even larger than for the Belle version (see Table HI]), and the r^ip') width, 

r( 7 V>') = 25 keV (set 6), r(jijj') = 60 keV (set 7), (48) 

is a bit small as compared to the estimates ([7]). The 7T + ir~J/ip data are described better 
than in the virtual-state version, while the description of the D°D*° data is rather poor (it 
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FIG. 4: Differential rates for the 7r + 7r~ J /tp channel (left plots) and D°D*° channel (right plots) 
(see Eqs. (|4"5j) and (JM|), respectively) with the parameters given by set 4 (upper plots) and set 5 
(lower plots). Parameters from these sets are presented in Table HH The distributions integrated 
over the bins, with resolution function taken into account, are shown as filled dots, experimental 



data (see Refs. [6 



34l |) are given as open dots with error bars. 



remains an open question either it is genuine, or is due to our ill-starred guess on the BaBar 
D°D*° resolution function). 

Therefore, we find that the BaBar data are more compatible with the assumption of the 
X(3872) being a virtual state of a dynamical nature. As to the charmonium admixture, we 
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FIG. 5: The same as in Fig. U]but for set 6 (upper plots) and set 7 (lower plots). 



can only state that it is small. Indeed, as shown in Ref. 41|, in the case of small values 
of W, the model independent Flatte analysis does not allow one to draw any conclusions 
on the binding mechanism, that is to distinguish between t-channel meson exchange forces 
or short-ranged s-channel forces due to coupling of bare states to the hadronic channel. 
One can only state that the properties of the resonance are given mainly by the hadronic 



continuum contribution, and the state is mostly of a dynamical 
We conclude this section with the comment on the paper 



molecular) nature. 



38 



where the Flatte fits 



were performed similar to ours, and the pole structure of the Flatte amplitude was studied. 
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Conclusions on the nature of the X were drawn in Ref. 
procedure developed in Ref. 



based on the pole counting 



46(: two near-threshold poles correspond to a large admixture 



of a bare state in the resonance wave function, while a single near-threshold pole indicates 
a dynamical nature of the resonance. Strictly speaking, the Riemann surface for the X case 
is much more complicated than the one assumed in Ref. [38], due to the presence of many- 
body cuts (caused by the n+n^J/ip and n + n^n J/ ; ip modes). However, the pole-counting 
procedure should yield, qualitatively, the same result as a more rigorous method based 
on the spectral density calculations employed here, which allows, inter alia, to estimate 
quantitatively the bare state admixture (for more details on the interrelation between the 
pole-counting and spectral density behaviour see Ref. (m])- For the fits with reasonably 
small values of the factor B two near-threshold poles were found in Ref. 38J, signalling a 
large admixture of the genuine charmonium, similarly to our Belle parameter sets. Note, 
however, that the fits presented in Ref. 38j are the overall ones: the Belle and BaBar data 
were fitted simultaneously. As a result, rather poor description of the BaBar D°D*° data 
was obtained, reflecting incompatibility of the new Belle and BaBar data. 



IV. COMMENT ON THE D* FINITE WIDTH 



The possibility of the bound-state solution brings on board one more important question. 
Namely, in the present analysis, the D*°-meson was assumed to be stable. As argued in 
Ref. 35], account for a small finite width of the D*° does not change the D°D*° line-shape 
in the case of the virtual state while, for a bound state, the effects of the finite width could 



be pronounced, as shown in Ref. 47( • A refined treatment of the finite width is in progress 
now [48], while here we estimate these effects using a simple ansatz suggested in Ref. 49j and 



re-invented in Ref. [47] . The recipe is to make the following replacement in the expressions 
for the D°D*° momentum entering the formulae for differential rates: 



Q(E)h(E) -> vW v/£ 2 + Tfji + E 



and 



G{-E) Kl (E) - v /^ V /v / £ 2 + ^74- E 
where i"! is the width of the D*°-meson. It can be shown 



(49) 



(50) 



48] that these formulae are valid 
if the resonance is well-separated from the three-body threshold (the D°D*°it threshold in 
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our case), and the zero-width limit is readily reproduced as r* — ► 0. 

To access the role of the finite D*° width we evaluate the D°D*° differential rates, with 
A = 63 keV, for the sets 2 and 3 (bound state scenario for the Belle data) and for the 
sets 4 and 5 (virtual state scenario for the BaBar data) and plot them, together with the 
zero-width rates, in Fig. El In particular, in this figure, we show our theoretical curves given 
by expressions (I23l) . with the replacement (j4"9|) . (|50|) . and with the Flatte parameters from 
the corresponding tables, without signal-background interference and not smeared with the 
resolution functions. 

As seen from Fig. El the virtual-state solutions are not affected by the finite D*° width at 
all (lower plots). In the meantime, the bound-state excitation curves are not affected either, 
in the above-threshold region, and a non-negligible D°D°7r° peak is developed around the 
bound-state mass position (upper plots). 

The bound-state peak resides at about —0.75 MeV for set 2 and at about —0.5 MeV 
for set 3, so the only effect expected is an increase of the number of events in the first 
near-threshold bin. Indeed, for the Belle bound-state solutions, we have calculated the ratio 
Ni/Ni of the number of events in the first (i = 1) and second (i = 2) non-empty Belle bins, 
with (iVj) and without (iVj) inclusion of the finite width: 

Nx/Ni = 4.31, N 2 /N 2 = 1.01 (set 2), 
Ni/Ni = 1.99, N 2 /N 2 = 1.00 (set 3), 

where the ratios above are calculated without resolution and signal-background interference. 
Clearly, a large value of the ratio N\/Ni does not cause problems, as Belle bound-state 
solutions underestimate the number of events in the lowest bin only (see Fig. [2]), and the 
number of events in higher bins is not affected by the finite-width effect. 

Thus a considerable number of the D°D*° events is to appear below the nominal D°D*° 
threshold in the bound-state case. In the meantime, the present experimental situation does 
not allow one to identify the bound-state peak. This must be attributed to the peculiarities 
of the data analysis: both BaBar and Belle Collaborations assume that the D°D ir events 
come from the D°D*°, distorting in such a way the kinematics of the below-threshold events 
and feeding artificially the above-threshold region at the expense of the below-threshold one. 



In a quite recent paper [50j the D°D°ti distributions were obtained with the recipe 



SJ), fl50|) . and, in order to describe the D°D*° data, the above-mentioned kinematical 
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FIG. 6: Upper plots: the differential rates for the D°D*° channel given by the Flatte formula with 
the finite width included (solid lines) and with the zero width (dashed lines) for parameters 
sets 2 (left) and 3 (right). Lower plots: the same as in the upper plots but for parameters sets 4 
and 5. 

distortion was corrected with some feedback from data processing arrangements. As a result, 
a nice description of the D°D*° Belle data was obtained with a bound-state solution (the 
description of the BaBar data is rather poor in Ref . [50] , quite similar to our BaBar bound- 
state solutions 6 and 7). Notice, however, that the solutions of Ref. \^ differ from ours 
in several respects. First, the data on 7i + 7i^J/ip and D°D*° modes are analysed separately 



in Ref. 



501 ]. so that it remains unclear whether the presented solutions provide a tolerable 
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overall fit. Second, there is no extra width Jq in the fits describing the D°D*° mode (in fact, 
there is no inelasticity at all in the best fits for the D°D*° data). Besides that, the scattering 
length approximation for the D°D*° scattering amplitude is employed in Ref. {50], which is 
not adequate for our solutions. 

Finally, related to the question of the finite D*° width is the problem of the interference 
in the decay chains X -► D°D*° -> D D°ti and X -> D°D*° -► D°D°ti . According to 
the estimates made in Ref. [40| . the interference effects could enhance the below-threshold 
D°D°n rate up to two times, however the effect is much more moderate above threshold 



48). As to the X -> D°D*° -> D°D°'y and X -> D°D*° -> D°D°'y decay chains, these 



are shown to interfere destructively 40fl. The proper account for the interference cannot be 



done in the over-simplified framework presented here, as this effect is to be included in the 



coupled-channel scheme from the very beginning 



48) 



V. DISCUSSION AND OUTLOOK 



The present analysis was prompted mainly by two recent experimental results: the dis- 

-iOn*o 



and the new Belle data |37J on the D°D 



36|. It is 



covery of the new r{^ijj') mode of the X(3872) 

mode. We have arrived at conclusions different from the ones of the papers [35j, 
instructive to discuss in detail the relation between these new results and the previous ones. 

To begin with, the r^jip') mode is to be included in the analysis. As the corresponding 
rate is comparable to the ir + 7T~J/if) one, virtual-state solutions yield a very large radiative 
decay width, about 500 -r- 800 keV which, at present, has no reasonable explanation. 

Furthermore, we have demonstrated that the new Belle data on the D°D*° mode are 
in conflict with the old Belle data [3J], as well as with the BaBar data on the same mode 
(up to the resolution issue, as was mentioned before). If the Belle peak at 3872.6 MeV 
is real, as suggested by the fine resolution and high statistics of the new Belle data, then 
there is no need anymore to reconcile the tt + it~J/iIj and DD* peaks, and the best fit is 
consistent with the bound-state solution. However, to overcome the problem of the large 
ratio (j34p of the branching fractions, we are forced to include "extra" non-DD* modes, with 
the radiative width r^ip') being only a small fraction of these "extra" modes. Thus the 



difference between the present results and the ones of Ref. [35J is due to the new data as 
well as the model-dependence allowed here. 
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The spectral density was calculated for all solutions presented and it appears that, for 
the bound-state solutions, there is a significant admixture of the bare state in the X wave 
function. The properties of the bare state delivered by the bound-state solution are in good 
agreement with the ones of the x' c t charmonium: the branching fraction in the B — > K decay, 
the total width, and the radiative jip' width comply well with the charmonium assignment. 

We stress that, in this picture, the X is not a bona fide charmonium accidentally residing 
at the D°D*° threshold. Had it been the case, the integral of the spectral density over the 
resonance region would have been unity while, for our bound-state solutions, it does not 



exceed 50%. It is rather a resonance attracted to the thresho 



in Ref. 43] and described in microscopical models in Refs. 21 



, a phenomenon advocated 



221 ] . In other words, the X 



is generated dynamically by a strong coupling of the bare x'ci state to the DD* hadronic 
channel, with a large admixture of the DD* molecular component. 

On the contrary, the virtual-state solution favoured by the BaBar data points to a rather 
small (if any) admixture of the bare state in the X wave function, and there is no need to 
invoke "extra" modes. This feature, in principle, could discriminate between bound-state 
and virtual-state solutions. In practice, the annihilation (light hadrons) modes encoded 
in the quantity r are not easily detectable so, in further studies, one is to rely upon 
improvements in the data on already observed modes. 

In particular, a clear signature for a bound-state solution is the below-threshold D°D 7t° 
peak. Unfortunately, from the experimental point of view, published data are not decisive, 
mainly due to the kinematical cuts imposed by the assumption on the D°(D )ir mode 
coming from D*°(D*°) one. In this regard, we urge both experimental collaborations to 
overcome this and to perform an unbiased analysis. 
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